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ABSTRACT 

Leading models of galaxy formation require large-scale energetic outflows to regulate the 
growth of distant galaxies and their central black holes. However, current observational sup- 
port for this hypothesis at high redshift is mostly limited to rare z > 2 radio galaxies. Here 
we present Gemini-North NIFS Intregral Field Unit (IFU) observations of the [O m]A5007 
emission from a z ps 2 ultraluminous infrared galaxy (Ltr > 10 12 L Q ) with an optically 
identified Active Galactic Nucleus (AGN). The spatial extent (« 4-8 kpc) of the high velocity 
and broad [O III] emission are consistent with that found in z > 2 radio galaxies, indicating 
the presence of a large-scale energetic outflow in a galaxy population potentially orders of 
magnitude more common than distant radio galaxies. The low radio luminosity of this system 
indicates that radio-bright jets are unlikely to be responsible for driving the outflow. However, 
the estimated energy input required to produce the large-scale outflow signatures (of order 
k 10 59 ergs over ra 30 Myrs) could be delivered by a wind radiatively driven by the AGN 
and/or supernovae winds from intense star formation. The energy injection required to drive 
the outflow is comparable to the estimated binding energy of the galaxy spheroid, suggesting 
that it can have a significant impact on the evolution of the galaxy. We argue that the out- 
flow observed in this system is likely to be comparatively typical of the high-redshift ULIRG 
population and discuss the implications of these observations for galaxy formation models. 
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1 INTRODUCTION 

The most successful current models of galaxy formation invoke 
large-scale energetic outflows to explain many of the properties 
of local massive galaxies and the intragalactic medium (IGM; i.e., 
red optical colours of massive galaxies, steep optical galaxy lumi- 
nosity function, the black-hole-spheroid mass relationship, metal- 
enrichment of the IGM; e.g., Silk & Rees 1998; Fabian 1999; King 
2003; Benson et al. 2003; Granato et al. 2004; Di Matteo et al. 
2005; Springel et al. 2005; Bower et al. 2006; Croton et al. 2006; 
Hopkins et al. 2006). A key attribute of these simulated outflows 
is the injection of significant amounts of kinetic energy into the in- 
terstellar medium (ISM), which can inhibit and terminate star for- 
mation by either heating the ISM or ejecting the gas out of the 
gravitational potential of the host galaxy. Large-scale outflows can 
be powered by star formation and/or Active Galactic Nuclei (AGN) 
activity (e.g., Heckman et al. 1990; Crenshaw et al. 2003; Veilleux 
et al. 2005), although most galaxy formation models predict that 



only AGN-driven outflows will be sufficiently energetic to have sig- 
nificant impact on the formation and evolution of massive galaxies. 
These large-scale outflows are expected to be particularly effective 
at z ~ 2, when star formation in the most massive galaxies was 
in significant decline (e.g., Juneau et al. 2005; Panter et al. 2007; 
Perez-Gonzalez et al. 2008; Damen et al. 2009). 

Broadly speaking, AGN-driven outflows are kinematically en- 
ergetic "winds" or "jets" that are initially launched close to the cen- 
tral black hole. The two main catalysts that have been proposed to 
power an AGN-driven outflow are (1) a radio jet or lobe, and (2) 
a radiatively driven wind. A radio-jet/lobe driven outflow will be 
most prevalent in radio-loud AGN but might only occur in a minor- 
ity of the AGN population^ For example, ~ 4-24% of the optically 
bright (radiatively strong) AGN population are radio loud, which 



1 The typical definition of a radio-loud AGN is vL5QH z /i^i440nm ~ 10 
(e.g., Kellerman et al. 1989). 
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appears to be a function of both redshift and AGN luminosity (e.g., 
Hooper et al. 1995; Jiang et al. 2007), although this may indicate 
the "duty cycle" of radio AGN activity. The fraction of optically 
faint (radiatively weak) radio-loud AGN fraction is a strong func- 
tion of stellar mass (e.g., Best et al. 2005; Smolcic et al. 2009). By 
comparison, a radiatively driven wind would be most effective in 
luminous AGNs with high mass-accretion rates, since the wind is 
directly driven by radiation pressure from the accretion disk. 

High signal-to-noise ratio X-ray and ultra-violet absorption- 
line spectroscopy have shown that high-velocity outflows are 
present in many AGNs, and may be a ubiqitious property of the 
AGN population at both low and high redshift (e.g., Crenshaw et al. 
1999: Chartas et al. 2002, 2007a,b; Laor & Brandt 2002; Pounds 
et al. 2003; Reeves et al. 2003; Porquet et al. 2004; Ganguly & 
Brotherton 2008; Gibson et al. 2009). For example, at least « 60% 
of unobscured AGNs show evidence for high-velocity outflows, 
and the maximum measured outflow velocity is a strong function of 
the AGN luminosity (e.g., Crenshaw et al. 1999; Ganguly & Broth- 
erton 2008). The estimated energy injection required to produce 
the X-ray absorption features (which sometimes suggest outflow 
velocities exceeding v ~ 0.1c) is typically ~ 0.1-1 of the AGN 
bolometric luminosity (e.g., Pounds et al. 2003; Reeves et al. 2003; 
Chartas et al. 2007a; Pounds & Reeves 2009). However, X-ray vari- 
ability studies and consideration of the required energy input to 
produce the high-velocity features, indicate that these outflow sig- 
natures must be produced in the vicinity of the accretion disc on 
< 1 pc scales (e.g., Crenshaw et al. 2003; King & Pounds 2003). 
Therefore, while it is clear that energetic outflows are present close 
to the accreting black hole of AGNs, it is far from clear what impact 
these outflows will have on the gas and star formation in the host 
galaxy on ~ 1-10 kpc scales. To directly test the impact that out- 
flows have on the formation and evolution of galaxies, it is neces- 
sary to identify large-scale energetic outflows, particularly in high- 
redshift (z ^ 2) massive galaxies where they are predicted to have 
been most prevalent. 

Optical spectroscopy of distant galaxies have revealed 
blueshifted absorption/emission lines in many systems, a signature 
of outflowing gas (e.g., Pettini et al. 2001; Shapley et al. 2003; 
Tremonti et al. 2007). However, the most direct insight into the 
identification and interpretation of large-scale outflows is provided 
by spatially resolved spectroscopy (integral-field unit; IFU obser- 
vations), which provide direct information about both the velocity 
and extent of any outflowing gas (e.g., Swinbank et al. 2005, 2006; 
Nesvadba et al. 2006). Using rest-frame optical IFU observations 
of distant radio-loud AGNs (i.e., high-redshift radio galaxies, here- 
after referred to as HzRGs), Nesvadba et al. (2006, 2007a, 2008) 
showed that large-scale energetic outflows are present in at least 
a fraction of the high-redshift galaxy population. The key diag- 
nostic that revealed an AGN-driven outflow in these systems was 
the presence of kinematically complex and extended [O III]A5007 
emission^ see Holt et al. (2008) for similar constraints on z 0.6 
radio galaxies. Nesvadba et al. (2008) showed that the extent of the 
broad [O III] emission is similar to that of the radio emission, pro- 
viding evidence for a causal connection between the radio emission 
and a large-scale outflow. The estimated kinetic energy required 

2 Emission from [O III] can be produced by a number of processes, in- 
cluding photo-ionisation and shocks (Osterbrock 1989) but it cannot be pro- 
duced in dense environments, such as the broad-line region of AGN, without 
being collisionally de-excited (e.g., Osterbrock 1989; Robson 1996). The 
detection of broad (FWHM > 500-1000 km -1 ) [O III] emission therefore 
indicates the presence of off-nuclear kinematic components. 



to drive the outflows in these HzRGs is m 1-40% of that poten- 
tially provided by the radio jet, indicating that they could be pow- 
ered by mechanical energy from the radio jet. However, HzRGs are 
rare systems (co-moving space densities at z « 2 of <E> pb 10 -8 - 
10 - 7 Mpc~ 3 ; e.g., Willott et al. 1998), possibly representing the 
most massive galaxies in the most biased regions of the distant 
Universe (e.g., Stevens et al. 2003; Seymour et al. 2007). To fully 
assess the global impact of large-scale energetic outflows on the 
formation and evolution of distant galaxies, and to distinguish be- 
tween the different mechanisms proposed to drive outflows, it is 
therefore necessary to search for large-scale energetic outflows in 
the more typical radio-quiet AGN population. 

In this paper we present Gemini-North NIFS IFU obser- 
vations of the [O III] emission from a z w 2 quasar (SMM 
J123716.01+620323.3, hereafter referred to as SMM J1237+6203) 
that is hosted in an Ultraluminous Infrared (IR) Galaxy (ULIRG; 
L m > 10 12 L Q ; e.g., Sanders & Mirabel 1996). ULIRGs (both at 
low and high redshift) are rapidly evolving galaxies undergoing in- 
tense dust-obscured star formation and AGN activity (e.g., Sanders 
& Mirabel 1996; Genzel et al. 1998; Page et al. 2001; Alexander 
et al. 2005b; Stevens et al. 2005; Schweitzer et al. 2006; Coppin 
et al. 2008; Lutz et al. 2008; Pope et al. 2008). Ultra-luminous 
infrared quasars are often believed to be a key phase in the evo- 
lution of a ULIRG, where energetic outflows are starting to shut 
down star formation in the host galaxy, before revealing an unob- 
scured quasar (e.g., Sanders et al. 1988; Canalizo & Stockton 2001; 
Granato et al. 2004; Hopkins et al. 2005a; Kawakatu et al. 2006). It 
is possible that every massive galaxy in the local Universe under- 
went a ULIRG phase at some time over the past ~ 13 Gyrs dur- 
ing which the galaxy spheroid and black hole were predominantly 
grown (e.g., Swinbank et al. 2006; Alexander et al. 2008; Tacconi 
et al. 2008). The [O III] luminosity of SMM J1237+6203 is compa- 
rable to the HzRGs studied by Nesvadba et al. (2006, 2007a, 2008), 
suggesting similar intrinsic AGN luminosities. However, impor- 
tantly, SMM J 1237+6203 is pa 3-4 orders of magnitude fainter at 
radio luminosities than the HzRGs. We identify the signatures of a 
large-scale energetic outflow in this ultraluminous infrared quasar, 
showing that these features are not restricted to the distant radio- 
loud AGN population. We estimate the energy input required to 
produce these features, assess what could be responsible for driv- 
ing this large-scale outflow, and discuss the implications for galaxy 
formation models. We have adopted Ho = 71 km s _1 , Q,m = 0.27 
and Q.a = 0.73; in this cosmology 1" corresponds to 8.5 kpc at 
z = 2.0. 

2 TARGET, OBSERVATIONS, AND DATA REDUCTION 
2.1 SMM J1237+6203 

The Gemini-North NIFS observations of SMM J 1237+6203 are 
part of an on-going NIFS program to identify large-scale energetic 
outflows in z ~ 2 galaxies bright at submillimetre (submm) wave- 
lengths0The optical position of SMM J 1237+6203 is q 20 oo = 12 h 
37 m 16=0, 5 2 ooo = +62°03'23'.'2 (Barger et al. 2003). 

SMM J1237+6203 has a submm flux of S 8 50=5.3±1.7 mjy 
(Chapman et al. 2005), is identified with a moderately bright VLA 
radio counterpart (Si.4=164.0±l 1 fjjy; Biggs & Ivison 2006), and 
has an optical spectroscopic redshift of z — 2.07 (Barger et al. 
2003). The rest-frame 1.4 GHz luminosity (fa 3 x 10 24 W Hz" 1 ) 

3 Programmes GN-2008A-Q-58 and GN-2009B-Q-21. 
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suggests a rest-frame far-IR luminosity of ~ 6 x 10 12 Lq, un- 
der the assumption of the local radio-far-IR luminosity correla- 
tion (Alexander et al. 2005a; e.g., Helou et al. 1985; Condon 
1992), identifying SMM J1237+6203 as a ULIRG. The detection 
of SMM J1237+6203 at submm wavelengths also suggests that it is 
bolometrically luminous (« 10 13 I/boi; e.g., Chapman et al. 2005). 
The strong detection of SMM J1237+6203 in the « l'.'5 VLA map 
(Biggs & Ivison 2006) but the non-detection in a sensitive (root- 
mean-square noise of 15.6 /iJy beam -1 ) C// 5 MERLIN map sug- 
gests that the radio emission is not dominated by a bright AGN 
core and is extended on scales of > 0'.'5 (R. Beswick, priv comm). 
This provides evidence that the radio luminosity (and therefore the 
far-IR luminosity) may be dominated by intense star-formation ac- 
tivity (~ 1000 Mq yr _1 ; calculated following Eqn. 4 in Kennicutt 
1998 using the radio-derived far-IR luminosity). 

The optical counterpart of SMM J 1237+6203 is bright (R = 
20.2) and optical-near-IR spectroscopy reveals the rest-frame UV- 
optical broad emission lines expected from a quasar (i.e., Lya; Nv; 
Civ; H<5-Ha with FWHM w 2100-2700 km s" 1 ; Chapman et al. 
2005; Takata et al. 2006; Coppin et al. 2008). Deep Chandra ob- 
servations reveal that SMM J1237+6203 is luminous (Z/2-iokcV ~ 



10 44 erg s 1 ) and possibly obscured (Nn 



")in the 



X-ray band (Alexander et al. 2005a). The apparently contradictory 
evidence of both luminous broad emission lines and X-ray absorp- 
tion imply that either the X-ray source and the broad-line gas are 
not co-spatial or that the X-ray absorption is due to an ionised out- 
flowing wind (e.g., Crenshaw et al. 2003); the latter would provide 
direct evidence for an AGN-driven outflow in the central regions, 
although higher signal-to-noise ratio data is required for confirma- 
tion. On the basis of the X-ray luminosity, the bolometric output of 
SMM J1237+6203 appears to be dominated by star-formation ac- 
tivity (Alexander et al. 2005a), which is consistent with the bright 
submm flux and the probable extended radio emission. However, 
we caution that we do not have unambiguous evidence for star- 
formation activity in SMM J 1237+6203; high signal-to-noise ratio 
spectroscopy would provide more direct evidence (e.g., the detec- 
tion of CO molecular gas or low-excitation star-formation emis- 
sion lines/broad-band features in mid-IR-mm spectroscopy). Low- 
resolution Subaru OHS spectroscopy shows that SMM J 1237+6203 
has luminous, broad, and extended [O III] emission (Takata et al. 
2006), providing evidence for the presence of a large-scale outflow. 

2.2 Gemini-North NIFS Observations 

The NIFS IFU uses an image slicer to take a 3'.'0 x 3'.'0 field with 
a pixel scale of 0'.'043 and divides it into 29 slices of width 0'.'103. 
The dispersed spectra from the slices are reformatted on the detec- 
tor to provide 2-dimensional spectro-imaging. For our observations 
we used the .ff-band grism at a spectral resolution of A/AA ~ 5290 
(a ^3A;v^ 25kms~ 1 ). They were taken on 2008 May 22 and 
2008 May 30 in excellent photometric conditions with ~ Of.' 3 see- 
ing. The observations were performed using the standard ABBA 
configuration in which we chopped by 6" to blank sky to achieve 
good sky subtraction. Individual exposures were 600 s and the in- 
tegrated exposure time was 15.6 ks (7.8 ks on-source and 7.8 ks 
on- sky). 

We reduced the data with the Gemini IRAF NIFS pipeline, fol- 
lowing the approach outlined in Swinbank et al. (2009). The data- 
reduction steps included sky subtraction, wavelength calibration, 
and flat fielding. We removed residual OH sky emission using the 
sky-subtraction routine described in Davies et al. (2007). To accu- 
rately align and and mosaic the individual datacubes, we created 
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Figure 1. Collapsed, one-dimensional integrated NIFS spectrum of 
SMM J1237+6203. The spectrum is well fitted by a two component model 
with both a broad and narrow component (dotted curves show the best fit- 
ting emission-line profiles), suggesting the presence of an energetic outflow. 
The inset plot shows the narrow [O III] velocity field extracted from the 
NIFS IFU datacube and plotted as a function of distance from the optically 
defined nucleus (region 4 in Fig. 2) along the major-axis. 



white light (wavelength collapsed) images around the redshifted 
[O III] line from each observing block and centroided the galaxy in 
the data cube. The data cubes were spatially aligned and co-added 
to create the final mosaic using a median with a 3a clipping thresh- 
old to remove remaining cosmetic defects and cosmic rays. Flux 
calibration was carried out by observing standard stars at similar 
airmass to the target galaxies, which were reduced in an identical 
manner to the targets. Since the [O III] emission lines lie in regions 
of > 99% sky transparancy, no corrections were required for tel- 
luric absorption. In all following sections, quoted line widths are 
deconvolved for the instrumental resolution. 



2.3 Emission-line Fitting 

We first constructed intensity and velocity maps of the [O III] emis- 
sion from the IFU data cube. We used a \ 2 minimisation procedure 
to fit each spectrum within the datacube, taking into account the 
greater noise at the positions of the sky lines. The spectra were av- 
eraged over a 3 x 3 spatial pixel (0'.'15 x 0'.'15), increasing to 5x5 
pixels and ultimately to 7x7 pixels if the signal was too low to 
give a sufficiently high x 2 improvement over a fit without the line. 
In regions where this averaging process still failed to give an ade- 
quate x 2 , no fit was made. To detect an emission line we required 
a A\ 2 > 25 (> 5 a) improvement in a continuum+emission-line 
fit over a simple continuum fit. When this criterion is met we fit the 
[O III] emission line with a single Gaussian profile, allowing the 
normalisation, width and central wavelength to vary. 

Since both broad and narrow [O III] emission-line components 
are seen in the collapsed one-dimensional spectrum (see Fig. 1), 
we also produced spatially binned spectra in 10 discrete regions 
to search for the presence of broad [O III] emission. For these in- 
tegrated spectra we fit both a single and double Gaussian profile 
emission line but only accepted the fit if the double Gaussian pro- 
file resulted in a significantly better fit than a single Gaussian profile 
(Ax 2 > 16 improvement; > 4 a). We present the emission-line 
fitting results in Table 1. 




Figure 2. [O III] intensity, narrow [O III] velocity maps, and line-free continuum image (centred at 1.51 fim; rest-frame 0.49 fim) of SMM J1237+6203 from 
the NIFS data cube. North is up and East is left. The circle at the top right-hand corner of the middle panel denotes the seeing disk for the observations and the 
grid of dashed lines indicate the sub-regions where spectra are extracted; see Fig. 3. The [O III]A5007 emission is evidently complex, with clear velocity and 
intensity gradients. The contours in the right-hand panel represent an intensity-weighted map of the broad [O III] emission components from Fig. 3, smoothed 
with a a = 2 pixel kernel (FWHM = 4.75 pixels), which has an angular extent of ss O'.'5-l'.'O (~ 4-8 kpc). 




Figure 3. Integrated spectra of SMMJ 1237+6203 for the ten sub-regions shown in Fig. 2. The red dashed curve denotes the best-fit emission-line parameters; 
in cases where a double Gaussian profile provide the best fit, both components are shown (dotted blue and magenta curves, offset in flux for clarity). The black 
dotted vertical line indicates the systemic redshift. The properties of the emission-line fits are given in Table 1 . 



3 RESULTS 

The collapsed one-dimensional integrated NIFS spectrum clearly 
shows strong [O III]AA4959,5007 emission; see Fig. 1. The 
[O III]A5007 emission line has a blue asymmetric profile, which is 
fitted with two underlying Gaussian profiles (the line morphology 
is also displayed in the weaker [O III] A4959 emission line). 
The narrow component of the [O III]A5007 emission is centered 
at a redshift of zn = 2.075 1±0.0001 (hereafter defined as the 
systemic redshift) with a width of FWHM = 212±28kms~ 1 . 
The broad [O III] emission-line component is centered at 
z B = 2.0725±0.0003 (Av = -254±60kms _1 from the narrow 
component) and has a width of FWHM = 823±118 km s _1 . 
Forty five percent of the [O III]A5007 emission is from 
the broad component, with an observed luminosity of 
^[oni] = (1-1 ±0.2) xlO 43 erg s _1 ; the observed narrow 



[O III] luminosity is L [0 iii] = (1.4±0.2) xlO 43 erg s . The 
broad [O III] luminosity of SMM J1237+6203 is comparable 
to the HzRGs studied by Nesvadba et al. (2006, 2007a, 2008; 
« [1-8] xlO 43 ergs -1 ). 

In Fig. 2 we show the [O III] intensity, narrow [O III] veloc- 
ity (measured with respect to the systemic redshift), and line-free 
optical continuum of SMM J1237+6203, measured from the NIFS 
data cube. The [O III] emission is strongly peaked and co-spatial 
with the optical continuum peak, as would be expected for an op- 
tically luminous AGN. However, there is also a bright spatially ex- 
tended narrow velocity component (~ l'.'7; ~ 14 kpc). The ve- 
locity map shows a strong gradient with a peak-to-peak shift of 
150±35kms _1 , with the highest and lowest velocity gas to the 
South-East and North- West of the Nucleus, respectively; see in- 
set plot in Fig. 1. If it is assumed that the velocity gradient of the 
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SMM J1237+6203 [O m]A5007 Emission-Line Constraints 



Region N comp Ax 2 2 (narrow) Flux (narrow) FWHM (narrow) 2 (broad) Flux (broad) FWHM (broad) A v (broad) 

(10~ 16 ergcm~ 2 s _1 ) (tar 1 ) (10~ 16 erg cm~ 2 s~ 1 ) (kmr 1 ) (kms -1 ) 
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0.41(6] 
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2 


45 


2.0751[5] 


0.25(7] 


132±33 


2.0740(12] 


0.58(5] 


468±85 
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50 


2.0750(4] 


0.51(7] 


209 ±42 


2.0728[(0] 


0.62(8] 


627 ±129 


224±96 
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2.0750(4] 


0.22(3] 


296±7( 










9 
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2.0752(5] 


0.28(2] 


275 ±66 











Table 1. The best-fit parameters to the [O III] A5007 emission line in SMMJ1237+6203. Values in [] denote the error in the last decimal place. The A% 2 
statistic corresponds to the improvement in \ 2 from fitting a double Gaussian emission-line profile over that of a single Gaussian emission-line profile. If the 
addition of the second Gaussian component gives a A% 2 > 16 (> 4 <r) improvement over the single Gaussian component fit, then we view the detection 
as significant and give the parameters of both (broad and narrow) emission lines. The regions (0-9) in column 1 denote those shown in Figs. 2 & 3. The Av 
velocity value corresponds to the offset between the systemic redshift and the wavelength of the broad [O III] component in each integrated spectrum. 



[O III] emission is due to host-galaxy rotation then the dynamical 
mass estimated from the velocity curve is fs 10 10 M for a canon- 
ical inclination of i — 30°, but could be as high as w 10 11 M© 
for an inclination angle of i = 10° (which is not potentially unrea- 
sonable given that this object is a quasar and therefore likely to be 
seen close to face on). However, due to potential complexities from 
dust extinction and broad [O III] components, it is not clear that the 
[O III] velocity field is dominated by galaxy rotation; see Lehnert 
et al. (2009) for other complications. By comparison, the spheroid 
mass estimated from the velocity dispersion of the narrow [O III] 
emission (« 200 km s" 1 ; see also Table 1) is of order « 10 11 M 
(calculated following Eqn. 1 of Erb et al. 2006 assuming C = 5), 
which is consistent with that expected from the estimated black- 
hole mass (M BH ~ 2 x 10 s M Q ; Alexander et al. 2008) given the 
local black-hole-spheroid mass relationship (e.g., Tremaine et al. 
2002; Marconi & Hunt 2003). 

In Fig. 3 we show the [O III] emission-line profile in the inte- 
grated spectra from the ten discrete regions (labelled 0-9 in Fig. 2). 
The [O III] emission is evidently kinematically complex in many of 
the regions across the galaxy. A broad (FWHM « 400-900 km s~ 1 ) 
[O III] component is identified in five (~ 50%) of the regions. The 
radial extent of the broad [O III] emission with respect to the nar- 
row [O III] emission is shown in Fig. 2 and Fig. 4. Four of the 
broad components (the exception is region 5, which has the lowest 
Ax 2 ) have a velocity offset (Av Si km s _1 ) with respect to the 
systemic redshift of SMM J1237+6203 and, interestingly, these re- 
gions are spatially coherent (over « 4-8 kpc) and are located to the 
South East of the nucleus; see Fig. 2. The presence of extended, 
high-velocity broad [O III] emission indicates a large-scale out- 
flow, due either to AGN and/or starburst activity. The lack of broad 
[O III] components to the North-West region of the nucleus could 
be due to extinction by the host galaxy, which may be expected if 
the outflow was driven by the AGN activity (i.e., the far side of the 
outflow would be obscured by the host galaxy). However, higher 
spatial resolution and signal-to-noise ratio data is required to pro- 
vide more definitive spatial constraints. 

4 DISCUSSION 

The Gemini-North NIFS observations of SMM J 1237+6203 show 
high-velocity (up-to « 300 km s _1 ), broad (up-to w 900 km s _1 ), 



large-scale (« 4—8 kpc) [O III] emission. The [O III] proper- 
ties of SMM J1237+6203 are consistent with those found for 
the HzRGs investigated by Nesvadba et al. (2006, 2007a, 2008) 
and indicate the presence of an energetic large-scale outflow in 
SMM J1237+6203. However, since SMM J1237+6203 is « 3^1 
orders of magnitude fainter at radio wavelengths than HzRGs, it is 
unlikely that the catalyst for this outflow is a radio jet. 

SMM J1237+6203 is the first high-redshift ULIRG with spa- 
tially extended broad [O III] emission to be mapped with IFU ob- 
servations. Previous IFU studies of high-redshift ULIRGs have typ- 
ically focused on the Ly a or H a emission line, and the majority 
of the objects observed have not shown evidence for AGN activity 
at optical wavelengths (e.g., Bower et al. 2004; Wilman et al. 2005; 
Swinbank et al. 2005, 2006; Nesvadba et al. 2007b). However, three 
pieces of indirect evidence suggest that SMM J 1237+6203 could be 
relatively typical of the high-redshift ULIRG population: (1) broad 
[O III] emission-line components have been identified with rest- 
frame optical spectroscopy in several high-redshift ULIRGs to date 
(e.g., Smail et al. 2003; Takata et al. 2006; Coppin et al. 2008), 
(2) w 50% of nearby ULIRGs hosting optical AGN activity have 
[O III] components with FWHM > 800 km s _1 , comparable to 
that found for SMM J1237+6203 (e.g., Veilleux et al. 1999; Zheng 
et al. 2002; see also Spoon & Holt 2009 for similar mid-IR spec- 
tral constraints), and (3) high-quality IFU data have been published 
for a number of nearby ULIRGs showing that they host broad and 
extended [O III], providing evidence for large-scale energetic out- 
flows in at least some ULIRGs in the local Universe (e.g., Colina 
et al. 1999; Wilman et al. 1999; Li'pari et al. 2009a,b). 

In this section we explore the properties of the large-scale out- 
flow in SMM J 1237+6203 in more detail by, first, estimating the 
energy input required to produce the observed outflow features (see 
§4.1) and, second, exploring what proceses could deliver sufficient 
energy to drive the outflow (see §4.2). Lastly, we consider the po- 
tential fate of the outflowing gas and briefly discuss the implica- 
tions of our observations for galaxy formation models (see §4.3). 



4.1 The Properties of the Large-Scale Outflow 

In Fig. 5 we plot the velocity offset and FWHM of the broad [O III] 
emission for SMM J 1237+6203 and compare it to HzRGs and sam- 
ples of 2 < 0.5 quasars. We have not corrected the widths and 
velocities of the [O III] emission-line gas for unknown projection 
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effects and dust extinction, and therefore they should be consid- 
ered lower limits. There appears to be a loose correlation between 
the width of the [O III] emission and the [O III] velocity offset, 
which would be expected for outflowing gas that becomes turbu- 
lent as the kinetic energy is dissipated through shocks into the ISM 
(e.g., Lehnert & Heckman 1996); however, we caution that this is 
unlikely to be the only possible interpretation. This figure shows 
that SMM J 1237+6203 has [O III] properties broadly similar to 
those of the HzRGs, with the exception of the most extreme source 
(PKS 1138-262) which has very broad [O III]; however, we have 
only plotted the integrated [O III] properties for the HzRGs here 
and note that the widths and velocity offsets of spatially distinct 
regions are often higher. This figure also shows that the integrated 
[O III] properties of SMM J1237+6203 are (1) not as extreme as 
some nearby ultraluminous infrared quasars, and (2) are more ex- 
treme than the majority of nearby quasars; indeed, only an extreme 
subset of the rare Narrow Line Seyfert 1 population (NLS1) are 
found to host [O III] properties similar to SMM J1237+6203 (e.g., 
on the basis of Komossa et al. 2008, only ~ 5-10% of NLS Is have 
high-velocity [Av £ 200 km s _1 ] broad [O III]). This further sug- 
gests that the origin of the outflow features in SMM J 1237+6203 
are coupled to the production of the infrared emission, which we 
explore in more detail in §4.2. However, we stress that the [O III] 
properties for the z < 0.5 quasar samples are from long-slit spec- 
troscopy and provide limited constraints on the spatial extent of 
these potential outflows. 



The extent of the broad [Om] emission-line gas from 
SMM J 1237+6203 is comparable to that found for the HzRGs stud- 
ied by Nesvadba et al. (2006, 2007a, 2008), which typically range 
from ~ 4-20 kpc, suggesting comparable overall energetics given 
the similar widths and velocities of the broad [O III] emission. 
Assuming that the broad [O III] properties are due to an energy- 
conserving bubble expanding into a uniform medium (see Eqn. 3 
in Nesvadba et al. 2006), we can place first-order constraints on the 
kinetic energy required to accelerate the broad [O III] emission- 
line gas to these velocities over large distances. For example, the 
estimated kinetic energy required to accelerate the broad [O III] 
emission to Av ~ 300 km s" 1 over ~ 4-8 kpc (the estimated 
size of the broad [O III] emission; see Fig. 2) is of order ~ (0.6- 
3) xlO 44 erg s" 1 for SMM J1237+6203, which should be com- 
pared to « (0.3-100) x 10 44 erg s" 1 for the HzRGs; in this calcu- 
lation we assumed the same ambient density as that estimated for 
the HzRGs by Nesvadba et al. (2008) since we do not have direct 
constraints for SMM J 1237+6203. 



Over a canonical 30 Myr quasar lifetime (e.g., Martini & 
Weinberg 2001; Hopkins et al. 2005b), the total injection of energy 
into the outflow would be of order w (0.3-3) x 10 59 ergs, which is 
comparable to the estimated binding energy of the galaxy spheroid 
in SMM J1237+6203; for example, based on an estimated spheroid 
mass of ~ 10 Mq and velocity dispersion of a ~ 200 km s _1 
(see §3), the estimated binding energy is ~ 10 59 ergs for R c = 
4 kpc (e.g., Binney & Tremaine 1987). This analysis is based on a 
simple model and should only be considered illustrative with uncer- 
tainties at the level of an order of magnitude but, given the limited 
constraints current available for high-redshift systems, a more com- 
plex model is not yet warranted. However, it does indicate that the 
large-scale outflow in SMM J 1237+6203 may be energetic enough 
to unbind at least a fraction of the ISM from the host galaxy. 



■ 400 

0) 

g 200 
o 
>> 
o 

■I o 
> 

o 

-200 

r-,1000 
I 

tn 

E 800 

S 600 

„ 400 h 
o 

— 200 



1 1 I 1 1 ' I ' 1 
- • Broad [OIII] 

o Narrow [OIII] 



I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 
SMMJ1237+6203. 



I I I I I I I I I I I I I I I I I I I 



i 4 



* * J J • * * 







-10 -8 -6 -4 -2 2 4 

R [kpc] 



Figure 4. Velocity and FWHM components of the broad and narrow [O III] 
emission of SMM J 1237+6203 plotted as a function of radius from the op- 
tically defined nucleus (region 4 in Fig. 2): the regions where broad [O III] 
emission is detected are labelled. The radius and bin width of the broad 
[O III] components are determined from Fig. 2 while the plotted narrow 
[O III] component is the narrow [O III] velocity field shown in Fig. 1. 



4.2 The Potential Power Sources of the Large-Scale Outflow 

On the basis of a comprehensive analysis of the broad [O III] prop- 
erties of HzRGs, Nesvadba et al. (2006, 2007a, 2008) argued that 
the outflows observed in HzRGs are most likely powered by ra- 
dio jets, which can potentially inject of order « 10 46 erg s _1 of 
mechanical energy into the ISM (although this calculation is de- 
pendent on the assumed radio luminosity to kinetic jet-power ra- 
tio). The kinematics, extent, and luminosities of the broad [O III] 
emission from SMM J1237+6203 are comparable to those of 
the HzRGs. However, since the 1.4 GHz radio luminosity of 
SMM J 1237+6203 is as 3-4 orders of magnitude lower than that 
found for the HzRGs, this strongly suggests that mechanical en- 
ergy from a radio jet is not the catalyst for the broad [O III] emis- 
sion seen in SMM J 1237+6203. For example, taking the same as- 
sumptions as those used in Nesvadba et al. (2006), we predict the 
mechanical energy from a radio jet in SMM J1237+6203 to be of 
order ~ 10 42 erg s _1 , which is two orders of magnitude below 
the required energy input to produce the observed broad [O III] 
features, even assuming a 100% energy-ISM coupling efficiency. 
This estimate of the radio-jet mechanical energy is also likely to 
be an upper limit since it is possible that a large, possibly domi- 
nant, fraction of the radio emission from SMM J 1237+6203 is due 
to star formation rather than AGN activity (see §2.1). However, 
could the high-velocity, broad, large-scale [O III] emission from 
SMM J 1237+6203 be powered by an outflow from a radiatively 
driven AGN wind or, given the indirect evidence for star-formation 
activity (see §2.1), supernova winds associated with an intense star- 
burst? 

Estimating the potential impact of an radiatively driven AGN 
wind on the ISM is difficult due to the huge difference in size scale 
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Figure 5. Broad [O III] emission-line velocity offset versus FWHM found 
for SMM J1237+6203, the HzRGs from Nesvadba et al. (2006, 2007a, 
2008), z < 0.3 ultra-luminous infrared quasars (Zheng et al. 2002), and 
samples of z < 0.5 quasars (Boronson et al. 2005; Komossa et al. 2008). 
The data for the nearby samples is from long-slit spectroscopy and provide 
limited constraints on the extent of the [O III] emission. The [O III] emis- 
sion from the HzRGs is extended on ^ 1 kpc scales but here we only plot 
the collapsed one-dimensional IFU spectra; however, for SMM J 1237+6203 
we show the data from the collapsed one-dimensional IFU spectrum (open 
circle) and the four individual regions that have broad [O III] emission with 
At; ^ km s — 1 (filled circles): the individual regions are labelled. The 
maximum velocities (i> ma x = v + 1/2 FWHM; i.e., Martin 2005; Rupke 
et al. 2005a,b) of the broad emission-line gas for SMM J1237+6203 are 
potentially high enough to exceed the escape velocity (v eBC ) of a massive 
galaxy (dotted curve), although since the majority of the gas has v < v csc 
it seems likely that only a minority of the gas will be expelled from the host 
galaxy. 



between the likely < 1 pc ejection radius of the wind and the large- 
scale (> 1 kpc) galactic enviroment. However, given the constraints 
from high signal-to-noise ratio X-ray spectroscopy for a number 
of AGNs to date (see §1), it seems likely that a significant frac- 
tion of the bolometric luminosity from the AGN will drive an out- 
flow close to the accretion disk, which at least provides a first-order 
constraint on the initial energy injection. Given the low signal-to- 
noise ratio of the X-ray spectra for SMM J1237+6203, we can- 
not directly identify the X-ray signatures of a high-velocity outflow 
in these data. However, we note that on the basis of high signal- 
to-noise ratio X-ray spectral analyses of the gravitationally lensed 
z = 3.9 ultraluminous infrared quasar APM 08279+5255 (poten- 
tially a higher-redshift analog of SMM J1237+6203), we may ex- 
pect » 0.1 of the AGN bolometric luminosity of SMM J 1237+6203 
to be injected into an outflow (e.g., Chartas et al. 2007a; see also 
Chartas et al. 2007b for similar X-ray spectral constraints on an- 
other distant ultraluminous infrared quasar); this energy injection 
fraction is also similar to that typically assumed in theoretical simu- 
lations (e.g., Wyithe & Loeb 2003; Di Matteo et al. 2005). The total 
bolometric luminosity of SMM J1237+6203 is ta 3 x 10 46 ergs -1 , 
although a significant fraction of this emission may be due to star- 



formation activity; see §2.1. A more conservative constraint on the 
AGN bolometric luminosity is ~ 3 x 10 45 erg s _1 , based on the 
rest-frame 2-10 keV luminosity and the average spectral energy 
distribution for quasars from Elvis et al. (1994). The initial energy 
input into the accretion-disk wind is therefore likely to be » (0.3- 
3) xlO 45 erg s _1 , which would correspond to a total energy in- 
jection over a 30 Myr quasar lifetime of ~ (0.3-3) xlO 60 ergs. 
Therefore, so long as ~ 10-100% of the energy from the wind can 
couple to the large-scale ISM gas, then the quasar can drive the out- 
flow to ~ 8 kpc (the maximum extent; see Fig. 2); the wind-ISM 
gas coupling efficiency drops to ~ 2-20% to drive the outflow to 
« 4 kpc (the minimum extent; see Fig. 2). These wind-gas cou- 
pling efficiencies are plausible given the likely large solid angle 
of radiatively driven AGN winds (e.g., Proga et al. 2000; King & 
Pounds 2003; Ohsuga et al. 2009). 

The dominant source of energy injection from intense 
star-formation activity is supernovae winds. Dalla Vecchia & 
Schaye (2008) estimate the energy input from supernovae to 
be « 10 49 erg Mg 1 , under the assumption of a Salpeter ini- 
tial mass function. Assuming a constant star-formation rate of 
~ 1000 Mq yr _1 (see §2.1), the energy injection into the host 
galaxy from supernovae winds would be » 3 x 10 44 erg s _1 . 
This suggests that the large-scale energetic outflow observed in 
SMM J1237+6203 could be powered by intense star-formation ac- 
tivity so long as ~ 20-100% of the injected supernovae energy 
can be coupled to the ISM (the range in wind-gas coupling frac- 
tions correspond to the range in [O III] extent of 4-8 kpc; see 
Fig. 2). Given that the energy injection from supernovae will be 
more widespread than that from the AGN (i.e., potentially over kpc 
scales rather than pc scales), such a large coupling fraction may 
not be unreasonable (e.g., Weaver et al. 1977). Furthermore, our 
calculation of the required energy to produce the broad [O III] fea- 
tures was based on the assumption that the outflow needed to be 
driven over ~ 4-8 kpc scales, when the distance could be substan- 
tially lower if the star-formation activity is widespread. For exam- 
ple, the energy injection required to produce the broad [O III] fea- 
tures over « 1-2 kpc scales is about an order of magnitude lower 
than that required to produce the broad [O III] features over ~ 4- 
8 kpc scales, which would give similar coupling fractions as those 
found in nearby starburst-dominated ULIRGs (e.g., Martin 2006). 
This suggests that even if the star formation in SMM J1237+6203 is 
lower than that assumed here then it is still possible to produce the 
broad [O III] features so long as the supernova winds can efficiently 
couple to the ISM. 

Both an AGN-powered and star-formation-powered wind 
are plausible candidates for driving the large-scale outflow in 
SMM J 1237+6203, and on the basis of our current constraints we 
cannot directly tell which process dominates. As mentioned above, 
broad [O III] components with FWHM > 800 km s" 1 are de- 
tected in w 50% of nearby ULIRGs hosting optical AGN activ- 
ity (e.g., Veilleux et al. 1999; Zheng et al. 2002) while, by com- 
parison, no nearby ULIRGs optically classified as HII galaxies or 
LINERs are found to host broad [O III] emission. This would seem 
to suggest that the production of broad [O III] is related to AGN 
activity. However, since [O III] emission is easily photoionised by 
AGN activity, the detection of broad [O III] components in opti- 
cally identified AGNs could be a sensitivity/selection effect (i.e., 
the AGN may just be more effective at illuminating the outflow- 
ing gas that is potentially present in all ULIRGs). Indeed, optical 
spectroscopy of the Na ID absorption line has shown that 50% 
of ULIRGs host outflowing neutral gas, irregardless of the optical 
spectral type (i.e., HII, LINER, Seyfert classification), indicating 
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that star-formation activity may often be responsible for gaseous 
outflows in ULIRGs (e.g., Martin 2005; Rupke et al. 2005a,b). It is 
therefore also interesting that typical quasars (i.e., those that are not 
infrared bright) rarely show such broad [O III] features (e.g., Boro- 
son et al. 2005; Komossa et al. 2008; see Fig. 5), which may reflect 
either an absence of intense star-formation activity or different evo- 
lutionary stages (i.e., an ultraluminous infrared quasar may be an 
earlier phase associated with energetic outflows; see §1). However, 
all of the evidence presented here is indirect and sensitive IFU ob- 
servations of larger samples of distant ULIRGs, AGNs, and star- 
forming galaxies are required to provide more direct constraints. 

4.3 The Potential Fate of the Outflowing Gas 

A crucial feature of large-scale energetic outflows in galaxy for- 
mation models is the potential to eject gas from the host galaxy. 
Although we have shown that the energy input from both AGN 
and star-formation activity is potentially large enough to drive the 
large-scale outflow and unbind the ISM, we have not investigated 
whether the gas can be ejected from the host galaxy. Murray et al. 
(2005) show that a significant fraction of the ISM can only be ex- 
pelled from the host galaxy if the energy injection exceeds a lim- 
iting luminosity, which is a function of the gas-mass fraction and 
the galaxy gravitational potential (this limiting luminosity is anal- 
ogous to the Eddington limit for the galaxy; see Eqn. 18 in Murray 
et al. 2005). We do not have an accurate constraint on the mass 
of the cold-molecular ISM gas in SMM J1237+6203; however, 
if SMM J 1237+6203 is similar to typical submillimetre-emitting 
galaxies then the moleculear gas mass would be of order « 10 10 - 
10 11 M© (e.g., Greve et al. 2005), which gives a gas-mass fraction 
~ 0.1-1.0. Assuming a = 200 km s _1 and a gas-mass fraction of 
0.1-1.0, the limiting luminosity for SMM J1237+6203 is « (0.8- 
7.8) xlO 13 Z/0, which is comparable to or larger than the total 
bolometric luminosity. Therefore, on the basis of this analysis, it 
seems unlikely that the majority of the ISM in SMM J 1237+6203 
will be expelled from the host galaxy. The velocities of the broad 
[O III] emission are also in broad agreement with this analysis, 
showing that while v mll , x >v esc in some of the regions, the typi- 
cal velocities of the gas are more modest (v < v esc ); see Fig. 5. In 
this respect, SMM J 1237+6203 is similar to nearby starburst galax- 
ies, which have maximum outflowing gas velocities within a tight 
range between the circular velocity and escape velocity of the host 
galaxy (e.g., Martin 2005). 

We therefore find that SMM J 1237+6203 is driving a large- 
scale energetic outflow that is disrupting the ISM and could po- 
tentially shut down star formation in the host galaxy, in qualitative 
agreement with predictions from models of the rapid growth phase 
of galaxies and their central black holes (e.g., Granato et al. 2004; 
Di Matteo et al. 2005; Hopkins et al. 2006). If SMM J1237+6203 
represents a typical ULIRG phase that all massive galaxies have 
undergone at some time over the past ~ 13 Gyrs, then this out- 
flow activity is potentially £s 2-3 orders of magnitude more com- 
mon than that found in HzRGs (i.e., the comoving space density of 
z ~ 2 submm-bright galaxies is $ w 10~ 5 Mpc~ 3 to be compared 
to $ « 10~ 8 -10~ 7 Mpc~ 3 for z « 2 HzRGs; Willott et al. 1998; 
Chapman et al. 2005; Swinbank et al. 2006). Indeed, we have iden- 
tified high signal-to-noise broad [O III] components in all of the 
three high-redshift ULIRGs observed in our on-going Gemini NIFS 
programs so far. It is therefore likely that the dominant outflow 
mechanism in distant galaxies is that found in SMM J 1237+6203 
(either AGN accretion-disk or supernovae winds) rather than radio- 
jet induced AGN activity. However, since the outflow observed in 



SMM J 1237+6203 requires significant gas accretion to drive either 
an AGN or supernova wind, once the quasar and star-formation 
activity has terminated the outflowing gas will cool, which will ini- 
tiate further star formation. Therefore, unless the gas entropy in 
the system has become so high that the gas-cooling timescales are 
> 1 Gyr (e.g., I. McCarthy et al. in prep), it seems likely that an- 
other process (e.g., a radiatively weak radio-loud AGN as is com- 
monly seen in nearby massive galaxies; Best et al. 2005, 2006) is 
required to keep the gas unbound and to prevent further star forma- 
tion from occuring on comparatively short timescales. 



5 CONCLUSIONS 

We have presented Gemini-North NIFS IFU observations of the 
[O III]A5007 emission in a z ~ 2 ultraluminous infrared galaxy 
hosting an optically identified AGN (SMM J1237+6203). Our main 
findings are the following: 

• The spatial extent (~ 4-8 kpc) of the high velocity (up-to 
v w 300 km s _1 ) and broad (up-to w 900 km s _1 ) [O III] emission 
is consistent with that found in z > 2 radio galaxies, suggesting the 
presence of a large-scale energetic outflow in a galaxy population 
potentially orders of magnitude more common than z > 2 radio 
galaxies. See §3, §4.1, & §4.3. 

• The estimated energy required to produce these large-scale 
outflow features (« [0.6-3] x 10 44 erg s _1 ) could be provided by 
a wind radiatively driven by the quasar (~ [0.3-3] x 10 45 erg s _1 ) 
and/or supernovae winds from intense star formation (~ 3 x 
10 44 erg s _1 ), so long as the wind-gas coupling efficiencies are 
comparatively high (of order £s 10-100%). However, the low radio 
luminosity of this system indicates that radio-bright jets are un- 
likely to be responsible for driving the outflow, in contrast to z > 2 
radio galaxies. See §4.1 and §4.2. 

• The estimated energy injection required to drive the large- 
scale outflow is comparable to the estimated binding energy of the 
galaxy spheroid in SMM J 1237+6203, suggesting that the outflow 
can have a significant impact on the evolution of the galaxy. How- 
ever, on the basis of the maximum energy input into the outflow and 
the measured outflow velocities, it seems unlikely that the majority 
of the ISM will be expelled from the host galaxy. This may mean 
that another process (or longer lived AGN activity) is required to 
keep the gas unbound and prevent future star formation from oc- 
curing. See §4.3. 

Finally, we conclude that SMM J1237+6203 is the first high- 
redshift ULIRG with spatially extended broad [O III] emission 
to be mapped with IFU observations, revealing the signatures ex- 
pected for a large-scale energetic outflow. However, since « 50% 
of nearby ULIRGs hosting optical AGN activity also have broad 
[O III] emission, we expect SMM J1237+6203 to be compara- 
tively typical of the high-redshift ULIRG population. High signal- 
to-noise ratio (and high spatial resolution) IFU observations of fur- 
ther high-redshift ULIRGs will provide more definitive constraints. 
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